N PS ARCHIVE 
1967 
BOYLE, R. 



STRUCTURAL DESIGN 


OF 


TRANSV 


ERSE 


RING 


FRAMES 








by 








Ron fid 


A. 


Boyl= 






Thesis Supervisor 






19 May 19^7 


J, Harvey Evans 











STRUCTURAL DESIGN OF TRANSVERSE RING FRAMES 



BY 



RONALD ANTHONY BOYLE 

it 

LIEUTENANT COMMANDER, UNITED STATES NAVY 
B.S., Marquette University 
(1958) 

Submitted in Partial Fulfillment of the 
Requirenents for the Degree of 
Naval Engineer and the Degree of 
Master of Science in Naval Architecture 
and Marine Engineering 
at the 
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
May, 1967 



Signature of Author: 



Department of Naval Architecture 
and Marine Engineering, May 18, 1967 



Certified by: 



Thesis Supervisor 



Accepted by: 



Chairman, Dernr rvntsl Connittee on 
Graduate Students 



LIBRARY 

NAVAL POSTGRADUATE SCHOOL 

MONTEREY, CALIF. 93940 



STRUCTURAL DESIGN OF TRANSVERSE RING FRAMES 



by 



LIEUTENANT CORIANDER RONALD ANTHOITY BOYLE 
UNITED STATES NAVY 



Submitted to the Department of Naval Architecture and Marine 
Engineering on 18 May 1967, in partial fulfillment of the 
requirements for the degree of Naval Engineer and the degree of 
Master of Science in Naval Architecture and Marine Engineer, 

ABSTRACT 

The relation of the total structure weight of a 
two-diir.ensianal ring frame to the stress distribution in each of 
the elements is investigated. The total weight of the two- 
dimensional structures investigated were a minirrun when design 
stresses had been attained in all elements of the structure. 

An expression is developed that approximates the change 
in stress distribution in a member when additions or deletions of 
material are made in the member itself and its associated members. 

This expression is used as the basis for a rational 
design technique that converges to a ring frame of minimum total 
weight. 
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I. INTRODUCTION 

The transverse structure of a wall sided shin consists 
primarily of bottom and deck beams, side frames and stanchions. 
This structure is called the transverse ring frame. The loads this 
structure is subjected to are side and bottom water pressure and the 
ship weights as well as dynamic loadings, 

The structural analysis of this type of structure, in two 
and three dimensions, is discussed in considerable detail in the 
literature. Computer solutions for both the two and three dimensional 
aspects are available, These rapid analysis tools are necessary pre- 
requisites for a rational design synthesis. 

In its most general form the optimization of the transverse 
structure entails the specification of the transverse frame spacing, 
the effect of longitudinal frame spacing, and the distribution of 
material throughout the frame. This thesis deals with the last part 
of the process, the optimium distribution of material throughout the 
ring frame of a wall sided ship once the transverse and longitudinal 
frame spacings are fixed. The optimum structure is defined as the 
least weight combination of deck beam, side frames, etc, that will 
carry the design loads. 
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II. MINIMUM WEIGHT STRUCTURES 

The end and field moments that a loaded bean element develops 
are directly related to its end restraint. These moment distributions 
are subject to wide variations. The bending moment diagrams for a 
beam element supporting a distributed load are reproduced in Figure 1, 

Full rotational restraint at both ends produces a symmetri- 
cal moment distribution and one in which the maximum moments occur 
at the ends, The field moment dominates when the moment restraint 
at the ends is zero. The intermediate restraiiit case is character- 
ized by a symmetrical moment distribution and identical end and 
field moments. In every case, limiting or design stresses are 
reached only at discrete points in the element. 

The size or weight of any element depends on the maximum 
moment developed when the design load is applied. The wide variation 
in maximum moments is accompanied by corresponding variations In the 
weight per unit length of members. Table 1 lists, for the ideal 
section developed in Appendix C, the variation in weight for the 
structural configurations of Figure 1, The absolute minimum weight 
per unit length occurs when the maximum stress occurs at three points 
simultaneously in the beam element, The specific weight increases 
by thirty percent going from the fixed-fixed beam to the free-free 
or fixed-free beam. Doubling the maximum moment for this ideal 
section increases the section weight by sixty percent. 
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The weight ratios of Table 1 underline the fact that the 
total weight of any structural system is related to how the individ- 
ual elements interact, A. minimum total weight of structure can be 
obtained by substituting members whose moduli satisfy the inter- 
mediate restraint condition. In most cases, the actual moment 
distributions that obtain when the structure is assembled will not 
be identical to those assumed conditions. The final total structural 
weight will depend on how the material distribution has been altered 
to obtain convergence. This, in turn, depends on the design criteria. 
The examples that follow illustrate how the total weight varies when 
the design criteria is that design stress must exist in at least one 
element; the stress in all other members must be equal to or less than 
the design stress. 

The first example, Figure 2, contains only elements sub- 
jected to bending. The length to depth ratio cf each element is 
assumed to be large enough to exclude shear effects. 

An initial choice of members based upon the intermediate 
restraint condition and whether the elements carry an external load 
gives: 

1/16 W ii 

Z 2 = -T— ~ l Z l = ° = Z 3 

Design 

The structure made up of these members does not possess adequate 
strength. The maximum stress in member two is greater than the 
design stress. The minimum weight members do not form an acceptable 
structure. 
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This structure can only be brought up to an acceptable 
level of strength by adding material. The weight of the final 
structure satisfying the design criteria will vary considerably de- 
pending on where the material is added, There are two directions 
available. One is to add material so that the. moments at the nodes 
are very small. This will cause the field moment in member two to 
predominate, The other members will be very sm.aH,, The second 
approach is to increase the size of all the members. This will 
bring about redundant moments at the nodes which are substantial, 

The weight curves of Figure 3 illustrate the trade off 
possibilities that exist when all the members are varied. Large 
increases in the specific weight of members one and three are 
accompanied by relatively small increases in member one as the dis- 
tribution factor approaches one, These effects are summarized in 
the total weight curve, These curves, for all the cases looked at> 
have definite minimums. When L = 5L.. = 5L , the curve shown in 
Figure 3, the minimum total weight occurred when all member sizes 
were equal and the node moments were maximum, At this minimum point, 
design stresses were obtained in all members. 

The following example is identical to the preceding except 
that all members have to support an external load. The loading is 
such that the fixed end moment of the middle span is twice that of 
the adjacent members; load symmetry is maintained. 
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As before, the structure made up of least weight members, 
corresponding to the intermediate restraint condition, is not 
satisfactory. All members are cverstressed. Sufficient strength 
can only be obtained by adding material. 

Total weight curves for this structure are shown in 
Figure 4, Each of these curves have minimum points which occur at 
k « ,5. At this point, all members are of equal size and the node 
moments govern. It is obvious then that the minimum total weight 
occurs when design stress levels exist in all the members. At all 
other values for the distribution factor, one or two members have 
too much material. 

Only bending stresses were considered in these examples, 
The substitution of a vertical member at each of the supports re- 
quires that direct stresses be accounted for (see Appendix B) ♦ The 
total weight curve for this structure also has a unique minimum. 
This minimum occurs when design stresses exist in all members. 

The iteration technique developed in this thesis conver- 
ges on a structure which, when the design load is applied, develops 
design stresses in each camber simultaneously. These stresses may 
develop at the ends of the members or at any point along the 
members. The structures which this design process converges to 
represents the optimum distribution of material for the given loads 
and the fixed frame spacing. The total weight of these structures 
is a minimum. Any other combination of members will result in a 
higher total weight of structure; the stress in some or all the 
members will be less than the design stress. 
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Ill, CONVERGENCE PROCEDURES 

The properties of a two dimensional transverse ring 
frame of a wall sided ship ray be considered similar to those of a 
plane frame. The members of this plane frame are subjected to bend- 
ing and direct stresses; side sway is not considered. In its most 
complex form, this space frame is made up of many bays and decks. 

Picking member sizes that are adequate, i.e. possess 
enough strength to support the applied loads, can be accomplished in 
a straight forward manner, St. Denis [1], presents a rational 
approach to this problem. The refinement of these initial* choices, 
i.e. cutting the strength of each member until it just meets the 
design stress, is discussed in this section, 

A two dimensional space frame of more than one deck is a 
statically indeterninent structure. The stress distribution result- 
ing from the design loads depend upon the size of each member of the 
frame. Changes in these stress distributions are expected if the 
size of one member is altered. The effects of changing the size of 
a member in a reasonably balanced structure will be greatest in the 
immediate vicinity of the member. It is assumed In the following 
that the stress or moment distribution in a member is affected only 
by changes in members that join the defining nodes of the member. 
This is the same as saying that each member of the frame is the 
same as member 2 of Figure 5. The stress distribution in member 2 
is affected by size changes in any of the members shown. 
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The refinement problem is this: given a frame with a fixed 
layout and subjected to specified design loads, what distribution 
of material will produce stress distributions such that in each 
member the maximum stress equals the design stress. This means 
that in each member that has a maximum stress less than design stress, 
the section is too heavy. Material must be taken away or the stress 
distribution altered so that the member carries more of the load. 
If the maximum stress in the member is greater than design stress, 
the member is undersized and the opposite holds. 

An approach to the problem that is suitable to machine 
applications is a selective procedure. After each stress analysis, 
the member in the worst condition, over or under design stress, 
is selected. The size of the member is altered based upon existing 
maximum moment. The structure with the corrected member undergoes 
another stress analyses. This cycle is repeated until all members 
are at the design stress level, 

A more appealing approach is one in which changes can be 
made in all the members before each stress analysis. The results 
of each stress analysis are used as the basis for changing unsatis- 
factory members. 

The implementation of this series approach requires that 
.interaction effects be estimated. Unless account is taken of how 
the stress distribution changes when the distribution of material is 
altered, the system will be subjected to undesirable oscillations. 
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The member shown in Figure 6 depicts a general member taken 
from the two dimensional plane frame. It carries a distributed 
load, W, and has non-symmetrical end restraint. The stress dis- 
tribution is not symmetrical; node one is assured to be overstressed . 

The information required is hoi; will this moment distribution 
change when any or all of the members involved change. At each, node 
three conditions have to be considered, an increase, decrease, or no 
change in the joint rigidity. 

An increase in the size of all the members adjacent to node 
two will provide additional rotational restraint. This is accompan- 
ied by an increase in the moment at node two and a decrease in the 
maximum moment. If the changes in the members adjacent to node two 
were sufficiently large, no change in member one would be required 
due to the altered stress distribution. 

Increasing the size of member one alone will cause the 
maximum moment to increase. If the size of the new member is picked 
based on the existing maximum moment, then it also will be over- 
stressed when the structure is reassembled. 

Table 2 is a summary of the conditions that have to be 
accounted for. Column three has the information needed to make a 
rational choice of new members. An expected increase in the maximum 
moment indicates that a larger than usual correction should be made 
in the already overstressed member. This would be accomplished by 
adding an increment to the existing moment, - The new member size would 
reflect this added increment. 
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Additional tables would be required for each of the 
following: a field moment governing, end moment identical and govern- 
ing, and non loaded elements. 

The Hardy Cross method of moment distribution offers a 
direct method for obtaining numerical estimates of changes in maxi- 
mum moments due to shifts in material distribution about nodes. The 
expression that is derived below actually embodies the tables dis- 
cussed in the last section. 

A first approximateion to the moment M of Figure 5 can 
be. obtained by carrying out the Hardy Cross method of moment dis- 
tribution through the first unlocking of all eight joints. It is 
not necessary to require that the joints three through eight be 
fixed , 

First release the joints adjacent to node two and 
distribute the unbalanced moments, U, The carry over from joints 
four, six and eight to node two is: 

CO = - 1/2 k^ U 6 - 1/2 k g2 U 8 - 1/2 k„ U^ 

This carry over is added to the unbalanced moment at node two, U„. 
The amount of this summation that carries over to node one after 
node two is unlocked is: 

1/2 k 21 [1/2 (k 62 U 6 + k 8 , U g + k^ u 4 > - U 2 ] 

The carry over from the joints adjacent to node one is 
found in the same way. It is: 

CO - - 1/2 k U - 1/2 k ?1 U ? - 1/2 k n U 3 
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"The moment of node one in member 1-2, M , can be 
expressed in the following approximate form: 



M 12 * FLM l-2 " k 12 U l + 1/2 k 12 (k 51 U 5 + k 31 U 3 + k 71 V 



1/2 k^ U 2 + 1/2 fc 12 k 21 U 2 + 1/4 k 21 (k 62 U 6 + 

k 82 U 8 + k A2 V 



1/4 k 21 k 12 (k 62 U 6 4- k 82 U g 4- k A2 U 4 ) + 1/4 k 12 k^ U x 

1/4 k^ 2 k 21 U 



(1) 

If M is the maximum moment in the member, the change in 
the moment as the adjoining members and member two itself changes is 
desired. The quantities FEM. . and U. depend only on the loads and 
the geometry and therefore are constants. The expression for M 
can be looked upon as a function of the distribution factors, k... 
which in turn depend directly on the section properties, 



H 12 " f (k 12 W 



The change in II is 



dM, = .TT— — dk , + rr — ' dk 

12 ' 3k i2 ij J 



or 



3f 3 f 

AM__ = ■—- Ak.„ + ~ Ak.. 

12 3k 12 ag. ij 
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Applying this to equation (1) results in: 



hi, U ' - U l Ak 12 + m (k 51 U 5 + k 31 U 3 + k 71 V Ak 12 



+ 1/2 k 12 «t 51 U 5 + 4k 31 U 3 + Ak ?1 U ? ) - 1/2 U 2 ik 



+ 1/2 U, (k 12 4k 21 + Ak 12 k 21 ) + 1/4 (k 62 U 6 + k 82 C + 



\ 2 V «< 21 



+ 1/4 k 21 (4k 62 U 6 + Ak g2 U g + 4k 42 U 4 ) + 1/4 ^ 

■ . < k 12 ik 21 + ik 12 k 21> 

(2) 

An adequate comparison of the orders of magnitude of the 
different parts of this expression requires knowledge of the joint 
unbalance. If it is assumed, as it was previously, that the loads 
are reasonably balanced throughout the structure, the structure of 
Figure 5 is again applicable. Applied to this structure equation 
(2) reduces to: 

AM 12 * ~ U l Ak 12 " 1/2 U 2 Ak 21 + 1/2 U 2 (k 12 Ak 21 + Ak 12 k 21 } 

+ 1/4 U x (k 12 Ak 21 + Ak 12 k 21 ) 

(3) 
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This is the sane as neglecting the unbalance at joints 
three through eight. The inclusion of these terras in a computer 
program is straightforward, but of course, would require more 
computational time. This final aporoximateion accounts for changes 
in any of the seven members involved, A more convenient form for 
equation (3) is; 



AM.. = [-U. (1 - 1/4 k .) + 1/2 k. . U.] Ak. . 



+ 1/2 [- U. (1 - k. .) + 1/2 U. k. .] Ak.. 



(A) 



Equation (4) gives a numerical value for the expected 



change in M if the adjoining members and member itself change size. 
This is an approximation as to how the moment distribution in member 
ij is affected when any or all the members change size. This is a 
valid approximation for both loaded and unloaded elements. There 
is no restriction to members undergoing only bending stress; it is 
equally useful for members with combination stresses. If different 
types of materials are used, their effects are included in the 
distribution factors. 

The use of this expression to find changes in the field 
moment is a little more involved. Equation (A) has to be evaluated 
for each end of the member in question if the moment distribution 
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is not symmetrical and added to the existing end moments, These 
new moments combined with the existing end shears give an approxi- 
mation to the change in field moment, 

A word about signs is appropriate. Clockwise moments 
are positive, counterclockwise moments are taken to be negative. 
The signs that result when AM. . is calculated cannot be interpreted 
literally, A positive value for AM. . may indicate an expected de-' 
crease in the moment, M , . , It is necessary to check the sign of 
the existing moment and then interpret the sign of AM , 
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TABLE 2 



Cumulative Effects of Member Variations 



(2) 



End 
w/ ^ 

Mom) 



End 
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Explanation of Table 

This table is drawn up for Member Two of Figure Five. It is 
assumed that the moment distribution is nonsymmetrical and that 
Member Two is overs tressed, Column three gives the expected change 
in the maximum moment based upon the cumulative effects of member 
changes summarized In columns one and two. 

Column 1: The sign indicates how the overall stiffness of 
the joint is expected to change based upon the 
stress conditions in the adjacent members. 
Column 2: The size of Member Two is expected to increase 
Column 3: The expected change in the maximum moment based 
upon the data of columns one and two. 
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IV. THE DESIGN SPIRAL 

The result of the previous section is a technique for 
estimating moment changes. These moment changes are brought about 
when the members at a node change their size. 

This is of use in the design process because it presents 
a method of estimating not only the magnitude of changes in the 
stress distribution (maximum moment) , but also the direction of 
the changes, It permits a more rapid convergence to the final 
structure because each member of the structure can be modified (if 
necessary) on each pass. And also because the members are picked 
based upon conditions that will prevail; not those of the last 
stress analysis, It is essentially a method of predicting the 
final moment distribution and picking the members based on this 
distribution. If the moment changed predicted were exact instead 
of approximations there wouldn't be any iterative procedure 
necessary, 

A desirable aspect of this technique, is that it is self 
modifying when applied repeatedly. The predicted moment changes 
are updated each time a new member is picked. In this way oscilla- 
tions are avoided and smooth convergence is assured, It is not 
necessary to perform a new structural .analysis until this self 
modifying aspect is completed. This occurs when there are no 
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significant- differences in the predicted moment changes from one 
cycle to the next. 

The rational design process developed in the sections that 
follow converges to the option distribution of material for a trans- 
verse ring frame. The general procedure followed is to use the 
results of a stress analysis as the basis for modifying the material 
distribution in such a manner that design stresses will be attained 
in each member, The results of Section II indicate that a structure 
which develops design stresses in each member simultaneously is 
the least weight structure. 

Equation (4) is used in this process to make corrections 
to the results of the stress analysis, A change in the scantlings 
of any member of the ring frame changes the stress distributions. 
As pointed out in Section III and summarized in Table It, the 
effect of a change in scantlings depends not only on the magnitude 
of the change, but also the existing moment distribution. Once the 
maximum moment in the member is identified, equation (4) gives 
direct information as to how this moment will change when the 
scantlings of the members of the ring frame are revised. 

This correction to the existing moment distributions 
provides the means of facilitating rapid convergence to the optimum 
structure. Changes in member scantlings are based on these pre- 
dicted moments. The objective of design stress in each member of 
the ring frame is obtained much more readily. 
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• The use of equation (A) in the design process requires 
the standard output of a stress analysis, end forces, plus end 
moments. The process begins with the results of the stress analysis 
which used the initial choices of members . 

The computer begins with member one, The governing 
moment is ascertained and the stress level is computed based on this 
maximum moment. If the stress level is within a certain prescribed 
range of the design stress, the member is not changed. This pro- 
vides a means for ending the iterations. If the stress level is 
not satisfactory, a new member size is calculated and stored in the 
matrix, K, This procedure is repeated for each member of '"the 
structure. The result of this cycle is the matrix K. The members 
that make up K represent the first cut at bringing about the desired 
distribution of material. These members are fictitious in the 
sense that they are not used in an actual stress analysis. The 
properties of these members are based upon the existing stress 
distributions. 

The computer now starts to revise the members of matrix K; 
it starts with the first loaded member. The stiffnesses, k..,of 
the members at the defining nodes of the member to be changed are 
computed using the data of matrix K, This data is used to compute 
AM,,, or AH., or both depending on the stress distribution, If the 
field moment governs AM, . is computed for each node; when an end 
moment governs AM., is calculated for the governing node. 
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The maximum morrent of the member to be changed is modified 
by adding the correction factor AM... The properties of a new 
section are calculated based on this modified moment. These 
properties are inserted in the proper slot of matrix K, 

The computer chooses the next loaded member and repeats 
the above. The unloaded members are revised after new properties 
have been picked for all the loaded members. This requires fewer 
steps because the end moments always govern. 

The result after all the members have been looked at is 

i 

a new- matrix of members, K , A stress analysis could be initiated 

at this point. The results of the examples indicate, however, that 

a smoother convergence is obtained if this process is repeated, i.e. 

n 

find the matrix K , 

A stress analysis is now undertaken using the members of 

ti 

the matrix K , The results of this stress analysis are used as 

the starting point for a nev; design cycle. 

Flow charts containing the details of the design spiral 
are contained in Appendix A, Illustrations of the design spiral 
in action are contained in Appendix B, 
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V. RESULTS AND CONCLUSIONS 

A rational design procedure has been developed which 
permits the rapid synthesis of the transverse ring frame of a 
wall sided ship. The resulting transverse structure represents 
the optimum distribution of material^ the minimum wieght combus- 
tion of deck beams, side frames, etc., within the constraints of 
fixed longitudinal and transverse frame spacing, 

Detailed applications of this procedure are carried out 
in Appendix E, Each step follows the subroutines developed for 
computer adaptation of this technique. Mo attempt has been made 
to adapt the procedure to hand calculations, 

This method converges rapidly. Large member changes were 
required in each of the examples to attain the design stress in each 
member. Example I required only one additional stress analysis 
after the initial choice of members were analyzed; Example II re- 
quired two additional moment distributions. 

The speed with which a satisfactory ring frame is obtained 
depends directly on the following: the limits set within which the 
member stress must fall and the use of weighting functions in con- 
junction with the predicted moment changes. 

The limits set for member stresses in Appendix B were 
liberal. This, of course, was conducive to a rapid convergence, 
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These limits are not fixed, however, and can be set at any desired 
level, It is possible to put new limits on the members stress with 
each problem. 

The predicted moment changes in the examples were used 
as given by equation (4) , The true approximation in equation (A) is 
in the magnitude of expected moment change. The direction of the 
change, plus or minus, is exact if the change is greater than a few 
percent. In all the cases investigated, the magnitude of the change 
was always low. Convergence would have been more rapid if the 
changes were multiplied by a weighting factor. This indicates that 
weighting factors, if necessary, can be developed as experience is 
gained with the basic program. These weighting factors would become 
part of the basic program. 

The results of Example II of Appendix B indicates that 
members whose continuity is not required to carry the applied loads 
may be eliminated. This is easily avoided by applying nominal 
loads to these members, e,g, their own weight. In an actual 
structure, this will seldom be a problem. 

The use of the ideal section of Appendix C was a necessity 
in the development stages of this program. The implementation of 
this design procedure, i.e. the actual programming, does not re- 
quire this ideal section. If a tape of actual sections is available, 
the program can use it and still work. The effect this would have 
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on the speed of convergence would be mixed. The actual members in 
most cases have larger section moduli than those the program would 
pick, ' The use of actual sections requires^ir, addition, the develop- 
ment of a trial ard error solution in the subroutine COIiSTRSSS , 
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APPENDIX A 



An algorithm covering the complete design process is 
developed, Flow cherts covering the details of the individual steps 
of the algorithm are also presented. 

Design Algorithm 



Start with member one compare the maximum 
and design stresses. If the maxirum stress is 
within the prescribed range, insert the pre- 
sent member properties in matrix K, If maxi- 
mum stress unsat,, pick a new irenber. Insert 
its properties in the matrix K, 



Repeat the above for each member. If all 
trcmbers have satisfactory stress levels the com- 
puter stops - the design is completed, If any 
or all members were changed, the iteration 
process is begun by returning to the loaded 
member number one. 



Required stiffness factors are calculated 
using elements in matrix K, The moment change 
is determined using this data. The existing 
maxirum moment plus this moment change is the 
new desirr. moment, A new rc-j.T is oic':ed usi: 
this design movent and the direct stress. 
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The new member is placed in matrix K, 
This process is repeated for all members 
even if stress levels are initially satis- 
factory. The matrix containing these new 

i 
members is now designated K , 



The data of matrix K is now used to 

calculate the required stiffness factors. 

The process is repeated. The matrix 

u 
formed is K , 



A stress analysis is carried out: using 
matrix K , The results of this analysis 
are the starting point for a new design 
cycle, 



The flow charts are made up in the form of subroutines. 
They cover all the steps presented in the design algorithm. 

The first flow chart is subroutine Moment. It uses as an 
input the data of the preceding stress analysis. This subroutine 
primarily finds the maximum moment and calculates the maximum 
stress. When field moments govern the member is given a special tag. 

Subroutine Change checks the maximum stress in each member. 
It stops the iteration process if all members are satisfactory. 
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Subroutine Stiffness uses the data of matrix K and the 
data of initial member choices to calculate the expected changes in 
stiffness, It also determines the moment change and adds this to 
the maximum moment, The net; member parameters are obtained by 
calling the appropriate subroutine, i.e. Pick cr Comstress, This 
subroutine assumes that a member numbering system is available which 
facilitates picking the required data from the matrix K, etc. 

Subroutines Pick and Comstress both choose new member 
properties, Comstress handles members that have significant direct 
stress as well as bending stresses. 

The calling sequence for these, subroutines would be: 

CALL SUBROUTINE MOMENT 

CALL SUBROUTINE CHANGE 

CALL SUBROUTINE STII'FNESS 

CALL STRESS - or any suitable stress analysis program. 
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SUBROUTINE MOMENT 



N-e Total # of Members 

Set Member Counter, n = 1 



'END 



Ji 






>/ 



K m; u ;Tv 






^/ M_ • FEM 

^ND * J 



ij 





I 






X <- 












My<- 


-M. . - V. . * X 
IJ ij 







^< "em ' -r >4 



Ij'?Iax,Field F 



END ^ ij 



M "S ^ M 

ij .lax ' * END 



J 



>/ F . . /A . . : 1 o n . V^. 

n. xj ij Design/ = 



a ) <: 

ij Max, Direct 



M..)„ F.. 

_ xA t Max _Jj_ 

" *Z . . A . . 

IJ ij 




n < n + 1 



KH^ 



Stop. Transfer 

Control to Main Program 
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SUBROUTINE CHANGE 




</oT) )r : 1.1 c n . X 
= < ij Max Design > 



r^/,90 n : a ) M \_JL 

= \ Design y }laxy 



Puts Present Member 
Properties in K 



n < n 4 1 



~c< n 



D^ 



Stop Print Out 
Matrix K, 
Problem Completed 



If Member Flagged for 
Combined Stress, Sub- 
routine Corrstress 
Called, Otherwise 
Subroutine Pick Called, 



n < n + 1 



^rn 



Stop. Transfer 
Control to Main 
Program 
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SUBROUTINE STIFFNESS 



N - 


g Number of Members 




n < 1 ; F < 1 



-Yes 



-< ^Is M 



ember Horizontal 



.*o~ 



I. I I. I. 

SUM <- --- + —-— i- + —L + -J~^~ 

L . L . . . L . L . . 

1 l+l J J ~ 1 



k. . <- 



SUM 



Ak. . <_ k. . - k. . 



Same as Left Side 
Except Different 
Members , 



Repeat the above for 
Node j 



Yes 



Does Field Moment Govern 



No 



AM. . <^ 

1J 



H <- 



aq. (--) 



' . . + AM . . 
1J 1J 



x < — w. ./v. . 
1J lj 



>L< M. . - V. . * X 




Finds M. . )., 

ij Max 

AM. .) M < eq. (4) 

lj Max 

\[ ) <. — ::. ,) w + am. , 

ij Max ij .xax ij 



Subroutine Corns tress 
Called if Member Flagged . 
Otherwise Subroutine Pick 



n < n + 1 

r 



4<EII> 



D^ 



r + i 



{E^U 



= Vl. 



Z>- aL 



Stop, ' : cr ; '.v 
Ce - ~ - 

Control to St:-;.;.; 
Analysis Pro^r-o, 
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SUBROUTINE PICK. 



Z.. <r 



Design 



D.. < (Z. ./.032) 1/3 



ij 



(,052) D' 



I.. < (.0154) D 

_2J 



Pveturn Control to 
Calling Program 



SU ^ROUTINE C(V [STRESS 



- F. Jo x .052 



M. .)„ x 12 

a "T7o32 
o 



- °/„ + 



J& 



A 27 



E> «? 




- D /o - 4 



/TT77 5 



27 



D <?' A + B 



A. . <=- 
Z.. <-- 
I. . <?■ 

3-J 



.052 D' 



.032 D" 



— .0154 



Transfer Control to 
Calling Subroutine 
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APPENDIX B 

ESWTPLES 

EXAMP LE I 

The seven narhercd frame shown on page 45 is used to 
illustrate the optimization technique developed in the body of 
the thesis, 

A . S^eJ-ecti-oj^oJ^ Ij^^^^ 

' The initial set of member properties is part of the data 
input to the cor^puter. The following were selected for this 
example. 

(a) Z = 240 in 3 ; I = 2000 in 4 

(b) Z 1 = Z 3 = 120 in 3 ; I = I = 820 in 4 

(c) Z, = Z c = Z^ = Z., = 120 in 3 ; I, = 820 in 4 

4 5 6 7 4 

These properties i;ere obtained using the ideal section of Appenii:: C 

B. _ Stress Analysis 

The Hardy Cross method of morrent distribution was used to 
obtain end tr.oments and reactions, The results of this first analy- 
sis are: 

M 12 " M 21 ' ! Vh.x = 68A KIP " 

N 17 = M 28 = M 2 ) Max = W5 KIP FI 

M 13 = V>lax = 95 K ' : F A " 18 ° KIP FT 



C, Subroutine Cha nge 

(1) The maximum .stress in each member is compared to 
the design stress = 40 KSI 

(2) The size of each unsatisfactory member is changed , 
The maximum stress in each member is: 

a.).. = 34 KSI 
1 Max 

o,) M = 50 KSI 

7 Max 

o_)., = 15 KSI + 10 KSI = 25 KSI - a combined stress member, 

3 Max 

The new properties for each member are: 

M ) 
2 Max 1 4 

z = -i-iiS* « 204 in ; I. - 1600 in 

Design 
Z = 148 in 3 ; I = 1100 in 4 

Z. = 67 in 3 ; A, = 3.3 in 2 ; I. » 415 in 4 

4 4 4 

The remaining member sizes follow by symmetry. These proper- 
ties make up the matrix K, The subroutine PICK would be ined 
to obtain the properties of cambers one through three; sub- 
routine COMSTRESS would be called to obtain the properties 
for members four through eight, 

D . Sub r outline J5 1. if f ne ss 

i 

(1) Calculates new stiffness factors, k.. 

ij 

,(2) Calculates the change in stiffness factors with 

respect to the. initial merrber choices. 



- 37 



D . Sub r o u t ine S t i f f ne s s (Cont'd) 

(3) Daternines AM. . and the new factitious maximum moment 

(4) Picks new members based one the new maximum moment. 
This subroutine begins with the first loaded member, member 1 
and calculates k 



k = .310 k = 



k = .184 k =-0 



Ak - + .126 Ak = 



These results are used in equation (4) to obtain AM-,, 

AM = + 50.4 K' ; 11 = + 495 + 50.4 = 545 KIP FT 

Subroutine Pick now comes up with 

1 3 f 4 

Z = 153 in { I ? = 1140 in 

These member properties are immediately inserted in the 

matrix K. This provides a constant updating of the data. 

The computer proceeds thru the structure doing the 

above for each member 

Member 2 

t i 

k 12 = * 448 k 21 * * 452 

k 12 = .450 k 2] = .450 

» » 

Ak = -.002 Ak = + .002 

» 



AM 1 _ 2 = 1.1K*; - M 12 = 635 KIP FT 



z' - 204 in 3 ; I_ «* 1600 in 4 



- 38 - 



D. Subroutine Stiffness (Cont'd) 
Member 3 



k 17 = .310 k = 



k 71 " ° 



Ak 17 « + .126 Ak = 

i i 

A^ = + 50.4 KIP FT M ? = 545 KIP FT 



1 3 * 4 

Z ? = 153 in ; I = 1140 in 



Member 4 



k 13 = .117 k = 

k 13 = ' 18A k 31 = ° 

t i 

Ak 13 - .067 Ak 31 - 

t t 

AM = 27 KIP FT M =73 KIP FT 

1 3 * 4 

Subroutine COMSTRESS 





3 




._/ j.h , 


3 


«J -U V-/ U-IL 


Me 


[rber 5 
. i i 








f 




k 15 


— 


.121 




k 51 = o 




Ak 15 


= 


.184 




Ak 51 - 




ik ; 5 


= 


-.063 




k 51 = o 



AM 5 = -25.2 KIP FT M 5 = 75 KiP FT 



i 3 i 2 ' i 
Z< = 55 in , A = 7.2 in ; I = 310 in 



Members 6 and 7 follow by symmetry. 
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» 

The matrix K at this point becomes the matrix K , 

Subroutine Stiffness repeats the above process and produces matrix 

ii 
K , The results of this procedure are summarized below 



s 



Member 


Z 


1 3 

(in ) 


1 2 
A (in) 


I 


(in 


1 




167 






1300 


2 




204 






1600 


4 




50 


6.9 




280 


5 




50 


6.9 




280 



The matrix K is now the matrix of actual member 
properties. A stress analysis of the structure made up of these 
members is performed'. The results of this analysis are: 

M. = 595 KIP FT a, = 42 KSI = 1.07 a 

1 1 o 

M_ = 630 KIP FT o = 40 KSI = a 

2 2 o 

M, = ,41 KIP FT a, « 36 KSI = ,9a 
4 4 o 

Stresses and moments for other members follow by symmetry. 

Subroutine Change now finds all the stress levels 
satisfactory. The members of matrix K represent the optimum 
structure. The design spiral is complete, 
EXAMPLE 2 

This example was chosen to illustrate the procedure for 
a nonsymmetrical system and one in which the field moment governs, 
at least initially. 
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The frame analyzed is shown on page 46 . The optimization 
process proceeds as illustrated in detail in example one, 
A* .. _ JvJrJr.hfll .- JJ.rll^.g.r,-? *£? e r t i e s 

Z 2 = 240 in 3 ; I = 2000 in* 

3 2 A 

All other members - Z = 40 in ; A = 6 in ; I - 210 in 

B. Stress Analysis 

M » 16S KIP FT 

U 2 =780 KIP FT 

M. = 168 KIP FT 
4 

M «* 84 KIP FT 
C . S u t j^r ou^ ^i T\p. .C { \j\nge 

a i = 50 KSI z = 50 in 3 ; I = 290 in 4 

o 2 - 39 KSI ; M p governing , 

Z = 240 in ; I » 200 in 
O, = 50 KSI Z * 

Combined Stress , , , 

o 5 = 38.5 KSI Z^ = 83 in; A^ * 9.5 in ; I - 560 in'' 

' 3 "' ' 4 

Z 5 ~ 37 in ; A = 5,6 in ; I = 190 in 

D, Subroutine Stiffness 



1 

Matrix K 










■ 


Member 


Z 


' 3 

(in ) 


A 


' 2 
(in ) 


l' (in 4 ) 


1 




56.1 







320 


2 




203 







1700 


4 




140 




13.5 


1050 


• 5 




23 




4.1 


100 
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The field moment governed In member 2 at this point. 
However once the second set of members were obtained, the end 
moments predominated. The results of the second pass is 

IT 

matrix K 



"3 "2 "4 

Member Z (in ) A (in ) I (in ) 



1 


56.1 


— — — 


320 


2 


203 




1650 


4 


190 


16.5 


1500 


5 


15 


3 


60 



E . S tress Ana ly s i s 



M 1 = 530 KIP FT 



M 2 •» 655 KIP FT 



M^ =114 KIP FT 



M- = 
F. Subroutine Chanre 



a, = 38,8 KSI = .97 
1 o 



a « 30 KSI - .75 
I o 



a, = 48 KSI = 1.25 
4 o 



o 5 = 
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These stress values require another pass. The trend is 
obvious, Members two, four and six are taking all the load. 
The second iteration with subroutine Stiffness produces the 
following matrices. 

i 

Matrix K - second pass 



Member 


' 3 
Z (in ) 


' 2 
A (in ) 


I (i 


1 


17 


— -— 


40 


2 


203 





1650 


4 


262 


21 


2000 


5 





____ 


__~_ 




Matrix K - second pass 

" 3 

Z (in J ) 

192 
262 



(m ) 



21 



I (in ) 

1550 
2000 



M "646 KIP FT 



M 4 *= 646 KIP FT 



o - 39.5 KSI = .98 a 
z o 

a, = 29.6 + 11.1 = 40.7 KSI 
4 



1.025 o 
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These stresses are within the specified limits, 

ii 

The design spiral is terminated. The members of matrix K , 

second pass, form the optimum structure. 
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UA 



© 



12 KIPS/FT 



_¥ ±. 



1AA 



24 KIPS/FT 



-^ w 



© 



All L's equal 



L = 20' 
12 KIP3/FT 



x. & 



V © 



2 © 



© 



2 © 8 ; 



© 



rrrrn 



STRUCTURE FOP, EXAMPLE B-l 



FIGURE B-l 
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M * . i > V 



© 



7 © 



© 



24 KIPS/FT 



X JL £ 



5 

VU 1 ' V V 



© 



A_1A "1 



© 



L - L, = L c - 10' 
3 4 5 



L, = L-. = 20 ' 

D / 



2 © 8 



© 



\\M \\ 



STRUCTURE FOR EXAMPLE B-2 



FIGURE B-2 
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APPENDIX C 
THE IDEAL SECTION 

The section moduli of available steel sections listed in 
the Manual of_ St e e l Co n sj. r u c^ti. on do not vary in a continuous 
manner. The use of the least weight sections compounds the 
problem. Each section is useful for a wide range of section moduli. 

These facts make it very difficult, if one is using the 
sections listed in the AISC Manual, to pick out trends. Also the 
fact that there are no mathematical relationships existing between 
the important properties of the sections compounds the problems 
encountered when using actual sections, 

A mathematical section, i.e. one in which there are 
specific relationships existing between the section parameters, 
avoids these problems. The section developed for this thesis is 
outlined below. It conforms well to minimum weight sections listed 
in the AISC Manual for larger section moduli. At the lower end of 
the scale there is a greater divergence from actual sections, 

Let t f = flange thickness 

t = web thickness 
w 

D = depth of section 
w f = flange width 

The relationships chosen for the section are; 

t, " 2 t 
f w 

D =50 t 
w 

w = 20 t 
f w 
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The expressions for the area, section modulus and moment 
of inertia are: 

A = .052 D 2 
Z = .032 D 3 
I « .0154 D 4 

A log-log plot of these expressions is made on page 49, 
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